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ABSTRACT 

Optical/near- infrared (optical/NIR; OIR) light from low-mass neutron star X-ray bi- 
naries (NSXBs) in outburst is traditionally thought to be thermal emission from the 
O ' accretion disc. Here we present a comprehensive collection of quasi-simultaneous OIR 

and X-ray data from 19 low-magnetic field NSXBs, including new observations of three 
^ . sources: 4U 0614-1-09, LMC X-2 and GX 349-f 2. The average radio-OIR spectrum for 

' NSXBs is a « -1-0.2 (where oc v") at least at high luminosities when the radio jet 

is detected. This is comparable to, but slightly more inverted than the a « 0.0 found 
for black hole X-ray binaries. The OIR spectra and relations between OIR and X-ray 
fluxes are compared to those expected if the OIR emission is dominated by thermal 
emission from an X-ray or viscously heated disc, or synchrotron emission from the 
inner regions of the jets. We find that thermal emission due to X-ray reprocessing can 
explain all the data except at high luminosities for some NSXBs, namely the atolls 
and millisecond X-ray pulsars (MSXPs). Optically thin synchrotron emission from the 
•/^ ' jets (with an observed OIR spectral index of athin < 0) dominate the NIR light above 

I Lx ~ 10^^ erg s^^ and the optical above Lx ~ 10'^^ erg s~^ in these systems. For 

' NSXB Z-sources, the OIR observations can be explained by X-ray reprocessing alone, 

although synchrotron emission may make a low level contribution to the NIR, and 
could dominate the OIR in one or two cases. 



o 
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1 INTRODUCTION 

In low-mass neutron star X-ray binary (NSXB) systems, 
matter is accreted onto a neutron star from a companion 
via a disc which radiates in the optical, ultraviolet and 
X-ray regimes. Some of the accreting matter and energy 
in X-ray binaries (XBs) can be released from the sys- 
tem through coUimated outflows (jets). Evidence for jets 
associated wit h NSX B s date back more than a decade 
JStewart et al. 19931: Bradshaw. Fomalont fc Geldzahlerl 



119991 : iFomalont! Geldzahler fc Bradshawl l200ll ) but not 
until recently has the e vidence emerged f rom a ny waveband 
other than the radio. ICallanan et al.l (|2002| ) found that 
infrared K-hand flaring in the 'Z-source' NSXB GX 17+2 
could not be explained by an X-ray driven wind or repro- 
cessed X-rays, but shared many properties with the radio 
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(i.e. jet) variability previously seen in the same source. It is 
worth noting that the near-infrar ed (NIR) counterpart o f 
GX 13+1 is also largely variable l|Charles fc Navloij [l99^ ) . 
Recently, iMigliari et al.l (2006') for the flrst time spectrally 
detected optically thin synchrotron emission from the jets 
of an 'atoir NSXB 4U 0614-1-09 in the mid-infrared. 

In the optical regime, many spectral and timing 
studies of NSXBs have established the presence of an 
accretion disc reprocessing X-ray pho t ons to optical 



wave l engths (e.g. iMcClinto ck et al. 1979; Lawren ce et al 



19831 : iKong et all [2000. : .McGowan et al.ii2003i : iHynes et al 
20061 ). In quiescence, the companion star can come 
to dominate the optical/NIR (OIR ) emission (e .g. 



Thorstensen. Charles fc Bowveij 1 19781 : IChevalier et al.l 
19891 : IShahbaz. Navlor fc Charled [l99i). as is the case in 



high- mass X-ray binaries (e.g. Ivan den Heuvel fc Heis j 
I1972D and low -mass X-ray binaries that accrete from a 
giant star. X-ray reprocessing is generally thought to 
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dominate the QIR emission of n on-quiescent NSXBs (see 
van Paradiis fc McClintockl 1 19951 for a review; see also 



Chen. Shrader fc Livid 1 19971 : ICharles fc Cod 12006^. In 



low-mass black hole candidate X-ray binaries (BHXBs), 
additional OIR emission mechanisms have been ob- 
serve d, in particul a r from the viscously heated disc 
(e.R iKuulkersI 1 19981 : iBrocksopp et alj I2OO1I: iHoman et al. 
200E ) and f r om c ompact jets (e.g . Han fc Hiellming 



1992 : iFended I2001I: ICorbel fc Fended [20021: iHvnes et al 



20031: IBrocksopp. Bandvopadhvav fc FendeT |2004 
Buxton fc Bailvnll2004 Russell et al. 2006). 



Recently, an anomalous transient NIR excess has been 
observed in a number of millisecond X-ray pulsar (MSXP) 
NSXBs at high luminosities, which is equivocal in nature. 
The source rnost studied is SAX J1808.4-3658, for which 
IWang et all l|200ll ) found a NIR flux almost one order of 
magnitude too bright to originate from X-ray heating. The 
NIR flux density was comparable to a radio detection of 0.8 
mjy (with a fla t 2.5-8.6 GHz spectrum) seen one w eek after 
the NIR excess. ICreenhiU. Giles fc Coutures I (|2006l ) also re- 
ported an 7-band excess in a different outburst of the same 
source, which they attributed to synchrotron emission. In 
addition, a variable /- and i?-band excess in XTE J0929-31 4 
seen on the same day as a radio detection (|Giles et al.ll2005l). 
and a transient NIR excess in X TE J1814-338 (iKrauss et al.l 
I2OO5I ) and ICR J00291+5934 ([Torres et al1l2007l ) were afl 
interpreted as synchrotron emission from the jets in the sys- 
tems. The NIR excess appears to be ubiquitously absent at 
lower luminosities. 

Steady, partially self-absorbed jets probably exist in 
low-magnetic fleld (B < 10~" G) NSXBs in hard X-ray 
states ( Migliari fc Fenderll2006l : [Massi ,2006). These include 
'atoll-type' sources in the 'island' state, 'Z-type' sources in 
the 'horizontal branch' and possibly the 'normal branchQ 
and tra nsients at low accretion rates (Lx ^O.ll/Edd) such as 
MSXPs. lMassH I 20061) shows on theoretical grounds that the 
existence of jets in NSXBs depends on the magnetic field and 
mass accretion rate, and the conditions required for jet ejec- 
tion are probably fulfilled for most of the NSXBs (Z-sources, 
atolls and MSXPs) with know n magn etic field strengths. In 
Russell et al. (2006; hereafter IPapeij [3). we predicted that 
emission from the jets could contribute or even dominate the 
OIR light of NSXBs at high X-ray luminosities, for sources 
in a hard X-ray state. This is expected if LoPTjot oc (a 
theoretical relation supp orted by radio-X-ray observations; 
iMighari fc Fended l2006l ) and Lopt.xr oc LP^^ (where XR is 
X-ray reprocessing in the disc) , as is pred icted and observed 
Jvan Paradiis fc McClintocklll993 : |Pap er 11. In the sample of 
iPapeii U there were too little data at high X-ray luminosities 
to test this prediction for NSXBs. 

In this paper we analyse quasi-simultaneous OIR and X- 
ray data and OIR spectra of a large, comprehensive sample 
of NSXBs in order to constrain the Lx-dependent dominat- 
ing OIR emission processes and to test the above prediction. 
Data are collected from transient NSXBs: atolls, MSXPs 
and Z-sources. We use the technique (|Papeij l3) of comparing 
the relations between the quasi-simultaneous OIR and X- 
ray fluxes with the relations predicted from models of three 



^ These states correspond to th e position of the source in th e 
X-ray colour— colour diagram; see lHasineer &: van der KIislll989l . 



OIR emission processes: a viscously heated disc, an X-ray 
reprocessing disc and synchrotron-emitting jets. We also in- 
spect the spectral index of the OIR continuum, which differs 
largely between emission from the disc and from the jets. If 
the radio/OIR jet behaviour is ubiquitous in NSXBs, its 
properties will constrain many parameters in these systems, 
in particular the power of the jets (if the jet radiative ef- 
ficiency can be constrained) and the wavelength-dependent 
jet contribution at a given X-ray luminosity. 



2 METHODOLOGY & RESULTS 

A wealth of OIR and X-ray data from atolls, Z-sources and 
MSXPs were collected from the literature and converted to 
intrinsic luminosities (or monochromatic l uminos ities) in the 
same manner as described in Section 2 of |Papeij l3. Much of 
the data we use here were obtained for IPaperl J, however we 



have sought additional quasi-simultaneous OIR and X-ray 
data (and OIR spect ral ene rgy distributions; SEDs, which 
were not collected for lPaperl lj) for the purposes of obtaining 
a comprehensive sample for this work (see Tables 1 and 2 
for the references of the new data). 

Estimates of the distance to each source, the optical 



HI abso rpt 
F lPapeil lf a 



are listed in Table 2 of iPaped lll and Table 1 of this paper. 
No data were included for which the OIR contribution of 
the companion star could be significant (i.e. quiescence for 
most sources) unless this contribution is known. For these 
data we subtracted the star light contribution. For the Z- 
sources, OIR data were only included when the OIR fiuxes 
were significantly higher than the lowest measured for each 
source in each waveband, to ensure minimal contamination 
from the companion. When two or more OIR wavebands 
were quasi-simultaneous, OIR SEDs were produced in order 
to spectrally constrain the origin of the OIR emission. In ad- 
dition to the data collected from the literature, we obtained 
OIR photometry of three NSXBs using two telescopes; the 
observations and data reduction are discussed in the follow- 
ing subsections. 

While quasi-simultaneity was achieved in most cases 
with use of the RXTE ASM X-ray daily monitoring, this 
is not possible for some historical OIR observations. In 
these cases we include the data only in the OIR SEDs. 
Similarly, for some sources OIR-X-ray quasi-simultaneity 
was achieved, but there was only one OIR waveband 
available so OIR colours and SEDs were not obtained 
(e.g. GX 17+2). Quasi-simultaneous OIR-X-ray data 
were available for five sources (LMC X-2, Cyg X-2, 
Cen X-4, 4U 0614+0 9 and Aql X-1) before t h e ad- 
vent of the RXTE (iKaluzienski. Holt fc SwankI Il980l: 



ICanizares. McClin tock fc Griridlavl jl980l: ICharles et all 
T980; Bonnct-Bid aud et al.l Il989l: iHasinger et al.1 Il990l : 
[van Paradiis et al.llT990l : iMachin et al.lll990l ). 



2.1 Observations with the Danish 1.54-m 
Telescope 

VRIZ-hwA imaging of two Z-sources, GX 349+2 and LMC 
X-2, were taken in July 2001 using the camera on the Danish 
Faint Object Spectrograph and Camera (DFOSC). Table 3 
lists the observations used in this work. De-biasing and fiat 
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Table 1. Properties and data collected for the 11 NSXBs not included in |Paperl l3 and IGR J00291+5934, whose parameters have been 
updated since Paper I (see Table 2 of IPapeffl for the properties of the remaining NSXBs). 



Source 


Type 


Distance 


Period 




McB 


Av 


A^H / 


qcs 


At 1 


Fluxes - 


= alternative 




/ kpc 


/ hours 


/ Mq 


1 Mq 




lO^lcm-^ 


(band, 


days 


data 


name 




(rcf) 


(ref) 


(ref) 


(ref) 


(ref) 


(ref) 


ref) 




refs 


(I) 


(11) 


(III) 


(IV) 


(V) 


(VI) 


(VII) 


(VIII) 


(IX) 


(X) 


(XI) 


IGR 


MSXP 


2.8+1.0 


2.457 


1.4 


0.039- 


2.5+0.3 


4.64+0.58 




1.0 


1, 36 


J00291+5934 




(1) 


(13) 


(16) 


0.160 (16) 


(1) 


(1) 








LMC X-2 


Z-source 


50+10 


8.16 


~ 1.4 


~ 1.2 


0.15 


0.91+0.07 




0.5 


20, 37, 


= 4U 0520-72 




(2, 3) 


(14) 


(17) 


(17) 


(20) 


(29) 






38 


XTE J0929-314 


MSXP 


> 5 (4) 


0.726 




~ 0.008 


0.42+0.10 


0.76+0.24 




1.0 


21, 37 


= INTREF 390 




(--8 + 3) 


(14) 




(4) 


(21) 


(30) 








Cir X-1 


Z-source 


9.2+1.4 


398 




4+1 


10.5+1.5 


19+3 




0.2 


37, 39, 


= BR Cir 




(5) 


(5) 




(18) 


(22) 


(22) 






40, 41 


XTE J1701-462 


Z-source 


8.5+8.0 








9+4 


9+5 




0.5 


37, 42 






(6) 








(23)* 


(6) 








GX 349+2 


Z-source 


9.25+0.75 








5+1 


7.7+1.0 




1.0 


7, 37, 


= Sco X-2 




(7) 








(7) 


(31) 






38 


XTE J1814-338 


MSXP 


8.0+1.6 


4.27 




~ 0.5 


0.71+0.10 


1.63+0.21 




1.0 


19, 37 








(.5) 




(19) 


(19) 


(19) 








GX 13+1 


Z-source 


7+1 








15.3+2.3 






1.0 


24 


= 4U 1811-17 




(9) 








(24) 










GX 17+2 


Z-source 


8.0+2.4 








12.5+1.5 


15+2 




1.0 


37, 43 


= 4U 1813-14 




(10) 








(25) 


(32) 








HETE 


MSXP 


5+1 


1.39 




0.016-0.07 


0.89+0.22 


1.6+0.4 




2.0 


26, 37, 


J1900.1-2455 




(11, 12) 


(15) 




(15) 


(26)* 


(26) 






44, 45 


XTE J2123-058 


Atoll 


18.4+2.7 


5.96 


1.3 


0.60 


0.37+0.15 


0.66+0.27 


0.77 (R, 




35, 46 


= LZ Aqr 




(5) 


(5) 


(14) 


(14) 


(27) 


(27)t 


34, 35) 






Cyg X-2 


Z-source 


13.4+2.0 


236.2 


1.78 


0.60 


1.24+0.22 


1.9+0.5 




1.0 


47, 48 


= V1341 Cyg 




(5) 


(5) 


(14) 


(14) 


(28) 


(33) 









Columns give: (I) source names; (II) X-ray classification; (III) distance estimate; (IV) orbital period of the system; (V) mass of the 
neutron star in solar units (assumed to be ~ 1.4M0 if unconstrained); (VI) mass of the companion star in solar units; (VII) and (VIII) 
interstellar reddening in V-band, and interstellar HI absorption column (*j4v and ^A^h are estimated here from the relation 
A^H — 1.79 X \{r cm A\T\ iPredehl Xlchmit^lliil); (IX) the companion star OIR luminosity contribution in quiescence; (X) The 
maximum time separation. At, between the OIR and X-ray obs ervations defined as qua si-simultaneous; (X I) R eferences for th e 
quasi-simultaneous OI R an d X-ray fluxes collec t ed. R efere nces: (1) [Torres et al] ll20q7^: f2)lBovd et al] 1I2OOOI); (3)lKovacsl (|2000|); (4) 



[jo^l 



- , ., ■ — - . . .) IKov ^___. , , - 

iGallowav et al.l Jiooj): f5)ljonker fc Nelemanj ||2004|): f6)lKennea et al.l ll2006h: ('7)IWachter fc MargonI lll996l): C8) ISt rohmav cTet alj 
||2003|): (9) iBandyopadhyav et al.l (Il999l); (10) 'Kuulkers et al.' (2002"); (11) 'Kawai fc Suzuki' (20051); (12) iGallowav et al.l ll2005bl): (13) 
IShaw et al.[|20 05): (14) Rittcr & Kolb (2003); (15) Kaarct et al. (2006); (16) Gallowav ct al. (2005a|); fl 7) ICrampton et al.llTgga ): (18) 
Johnston, Fender & Wu (1999); (19) Krauss ct al. (2005); (20) Bonnct-Bidaud ct al. (1989); (21) Giles et al.' (2005"); (22) 



iJonkcr. Nclcmans h Bassa 12007 ); (23) Produit ct al. (2006 ); (24 ) Cha rles & Naylor (1992); (25) Dcutsch ct al. (1999); ( 26) 



ISteeghs et ah ll2005bl); (27)lHynes et al.l llioOll); (28) iMcClintock et al.l jl984i); (29) ISchu lz (1999 ); (30 ) |juet t. Gallowav & Chak rabart^ 
ll2003t); (31)llaria et atl ll2004i); (32) IVrtilek et aj] (Il99ll ); (33) ICostantini. Frevberg fcPredehli (120051); (34) iCasarcs ct al.. (.2002); (35) 



IShahbaz et al.l (|2003|); (36) ISteegshs et_ al.l (1200^(37) RXTE ^ ASM ; (38) this paper ; (39) IShirel (Il998l); (40) iGlassI ( 11973): (41) iMonetil 

al.l y 



( I1992I) ; (42) iMaitra fc Bailvnr(l20od) ; (4 3) ICallanan et al.l (|2002|); (44) (|2(305|); (4 5) 'Steeg hs et al.l (l2005a^; (46) iTomsick et al.l 
(|2004[); (47) lHasinger et al.ll990l); (48)lvan Paradiis et al.l (Il990l); ( 49) Tvan Paradiis ct al. (198d); (50 ) ICorbel et al.l (Il998l); (51) 
IWachtej (Il997[) : (52) iGilfanov et al.Hl998l ): (53) ICampana et al.l (|20M ): (54) iGreenhTil et al.l (|20o3 ): (55) iBailvn. Neil fc Maitral (|2006l ) 



fielding was performed with IRAF and aperture photometry 
of the targets and the standard star L TT 7987 (V = 12.23 
mag; R = 12.29; / = 12.37; Z = 12.50; iHamuv et aLlll992l : 
I1994I ) was achieved using PHOTOM. GX 349+2 was detected 
with a significance of > 80" in all exposures (most were > 
50cr detections). For LMC X-2 we aligned and combined 
the six exposures in each filter on MJD 52115 to achieve a 
higher signal-to- noise ratio (S/N). No combined images were 



created from the observations of this source on MJD 52117 
as the S/N was sufficiently high for photometry. LMC X-2 
was detected at a significance level of > 4(t in the Z-band 
images and > 20a" in all other images used in this work. In 
Fig. 1 we present the finder charts for LMC X-2. 

The measured fiuxes were accounted f or airmass- 
depeii dent atmospheric extinction according to iBurki et al.l 
(|l995l ). We fiux-calibrated the data using LTT 7987. For 
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Table 2. References of the new data from sources in |Papeil l3. 



Source 


Classifi- 
cation 


At 1 
days 


Fluxes - new 
data references 


IGR J00291+5934 


MSXP 


1.0 


1 


4U 0614+09 


Atoll 


0.5 


38 


CXOU 132619.7- 


unknown 






472910.8 








Cen X-4 


Atoll 


0.5 


49 


4U 1608-52 


Atoll 


1.0 


50, 51 


Sco X-1 


Z-source 


1.0 




SAX J1808.4-3658 


MSXP 


1.0 


37, 52, 53, 54 


Aql X-1 


Atoll 


0.5 


37, 55 



References: see caption of Table 1. 



GX 349+2 we measured the flux of three stars in the field of 
view with known V and it^-ban d magnitudes (stars 2, 5 and 
7 listed in lWachter fc Margon| [l996) and found their mag- 
nitudes to differ from those previously reported by ~ 0.04 
mag. The resulting fiuxes obtained for GX 349+2 and LMC 
X-2 were de-reddened to account for interstellar extinction 
using the values of A\i listed in Table 1 (Ar, A i and At, 
were calculated according to the recipe in iPaperl H. 

2.2 UKIRT Observations of 4U 0614+09 

NIR imaging of 4U 0614+09 was obtained with the 3.8 
m United Kingdom Infrared Telescope (UKIRT) on 2002 
February 14 ( MJD 52319.3), usin g UFTI, the UKIRT Fast 
Track Imager l|Roche et al.ll2003l ). Jittered observations of 
4U 0614+09 were made in the J, H and K filters, with 
9x30 s exposures in both J and //, and 9x60 s exposures in 
K. The infrared standard star FS 120 (J = 11.335 mag; H = 
10.852; K = 10.612) was also observed for a total of 50 s, 25 
s and 25 s in J, 7? and K, respectively. The airmass of the 
standard and 4U 0614+09 were very similar: 1.004-1.036. 
The 'JITTER-SELF-FLAT' data reduction recipe was used, 
which created a fiat field from the sequence of 9 jittered ob- 
ject frames and a dark frame. After dark subtraction and fiat 
fielding, a mosaic was generated from the 9 object frames. 

Photometry was carried out using IRAF. 4U 0614+09 
was detected with a significance of 7.1a in J, 15.3a in H 
and 11.5a in K. Flux calibration was achieved using FS 
120, yielding the following de-reddened {Ay = 1.41) flux 
densities for 4U 0614+09: = 0.145 ±0.037 mjy; F^,h = 
0.139 ±0.020 mJy; F^,k = 0.111 ±0.022 mJy (the apparent 
reddened magnitudes are J = 18.12; H — 17.50; K — 16.38). 

2.3 Results 

In the upper panels of Figs. 2 and 3 we plot the quasi- 
simultaneous OIR and X-ray data. The OIR monochromatic 
luminosity, L,/,oir is plotted against L^^a (where a is the 
orbital separation) in Fig. 2 in order to test the X-ray re- 
processing model. In Fig. 3, I/^.om is plotted against Lx: 
the plot necessary to test the models of jet and viscous disc 
OIR emission (Pap er I) . Thes e plots are the same as the left 
and right panels of Fig. 5 in IPaperl li with the new NSXB 
data added and the data from BHXB s remov ed. 

The orbital separations, as with iPaped U are inferred 
from the best known estimates of the orbital period and 




Figure 1. High resolution optical finding chart for LMC X— 2 
in y-band (upper panel) and J-ba nd (lower panel) (for a lower 
resolution S-band finding chart see Ijohnston et al.lll979l) . North 
is to the top and east is to the left. The images are ~ 0.8 X 0.8 
arcmins and were taken on 2001-07-27 with DFOSC on the Danish 
1.54-m Telescope. 

masses of the neutron star a nd companion (listed in Ta- 
ble 1, and Table 2 of |Papeij l3). For systems with obser- 
vationally unconstrained neutron star masses, we assume 
Mco ~ 1.4Mq. We do not include data in Fig. 2 from sys- 
tems in which the orbital period or the companion mass is 
unconstrained. We note that data from Cir X-1 in Fig. 2 
(open circles in the top right corner of the figure) may not 
be representative because its orbit is eccentric and so the 
orbital separation cannot be accurately inferred using this 
method. 

The power law lines in Figs. 2 and 3 repre- 
sent the expected correlations for X-ray reprocessing 
(|van Paradiis & McClintock' 'l994 }) an d jet emission (solid 
line; Migliari & Fender 2006; P apeil ||), respectively. The 
j et model is normalised to the observed radio-X-ray data 
(|Mighari fc Fender! I2OO6I ) assuming a fiat (optically thick) 
radio-OIR jet spectrum and the X-ray reprocessing m odel is 
normalised to the the optical (BVRI) NSXB data of iPapeJ 
d. The dotted power laws in Fig. 2 represent the expected 
relation for NIR data dominated by X-ray reprocessing as- 
suming an OIR spectral index of a = 0.5 and a = 2.0 (where 
Li, oc I/"). Although the different orbital inclinations be- 
tween these sources could change the X-ray and OIR lumi- 
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Table 3. Log of the Danish 1.54-m Telescope observations. 



MJD 




Target 


Exposures 


Integration time / exp. 


Apparent magnitudes (not de-reddened) 








/filter 


V 


R 


I 


Z 


V 


R 


I 


Z 


52114.07- 


-.25 


GX 349+2 


9 


300 


200 


120 


120 


18.46-18.62 


17.44-17.62 


16.54-16.65 


16.04-16.24 


52114.99- 


-5.07 


GX 349+2 


4 


300 


200 


120 


120 


18.32-18.57 


17.35-17.53 


16.47-16.58 


16.07-16.17 


52115.30- 


-.31 


LTT 7987 


1 


2 


2 


4 


6 










52115.35- 


-.43 


LMC X-2 


6 


120 


120 


180 


240 


19.19±0.07 


19.01±0.07 


18.91+0.07 


19.12+0.14 


52117.37- 


-.44 


LMC X-2 


2 


120 


120 


180 


300 


18.47±0.07 


18.43±0.07 


18.45+0.08 


18.57+0.26 



MJD 52114.0 corresponds to 2001-07-24.0 UT. The filters used were Bessel V, Bessel R, Gunn i and Gunn z. LTT 7987 is the standard 
star used for flux calibration. For GX 349+2 the range of magnitudes measured are tabulated, whereas for LMC X— 2 we tabulate the 

magnitudes from the combined images. 



nosities o bserved and hence these X-ray-OIR relations, we 
showed in |Papeij l3 that inclination doesn't appear to play a 
significant role. 

The lower panels of Figs. 2 and 3 describe the depen- 
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dence of the shape of the OIR spectrum with a and 
I/x, respectively. The spectral index a is calculated for all 
OIR data where two or more optical or NIR data points are 
quasi-simultaneous. The spectral index of the optical (us- 
ing data in the B, V or i?-bands) and NIR {R, I, J, H 
or Ti'-bands) data are shown separately. We define the R- 
band as the break between optical and NIR here because 
it is often the R — I colours that indicate the NIR excess 
(and for many of these observations JHK data were not 
obtained). For consistency (and to separate the two wave- 
bands used to calculate a) we use the two lowest frequency 
NIR bandpasses in each OIR SED to calculate cvnir and 
the two highest frequency optical bands to calculate aoPT. 
Data from the Z-sources are shown as different symbols to 
the atoUs/MSXPs because the spectral index of the former 
group appears to behave differently to that of the latter (see 
Section 3.2). 

We have also compiled luminosity-calibrated OIR SEDs 
of 17 NSXBs (Fig. 4). The references for the data are hsted 
in Tables 1 and 2, and Table 2 of Paper I. In Section 3 we 
attempt to interpret the OIR-X-ray plots and OIR SEDs in 
terms of the dominating optical and NIR emission processes. 



3 INTERPRETATION & DISCUSSION 
3.1 OIR - X-ray Relations 

In the upp er panels of Figs. 2 and 3 we improve on Fig. 5 
of iPaperl nl with a larger sample of NSXB data. Focusing on 
the upper panel of Fig. 3, we se e that the global relation 
Lorn oc L^^ for NSXBs (|Papei| [l does not hold at high 
X-ray luminosity (Lx > 10'^^'^ erg s~^); the OIR is more 
luminous than expected from this relation. This is inconsis- 
tent wit h the model for emission from a viscously heated disc 
(|Papeil llh. which requires a relation Loir oc for optical 
data and Loir Lx"^ for NIR (Fig. 3). We can therefore 
rule out a viscous disc origin to the OIR emission, at least 
at Lx > lO^''-^ erg s~\ 

The upper panel in Fig. 3 shows that the jet model lies 
close to the data at high Lx, suggesting the OIR enhance- 
ment at these luminosities may be caused by the domination 
of the jet. The jet model cannot describe the data at low Lx 
(< 10^^ erg s-i). However, the upper panel of Fig. 2 indi- 



cates that all the optical data also lie close to the expected 
relation for X-ray reprocessing in the disc. Hence, only the 
data above Lx > 10^^ erg s~^ could arise due to the presence 
of a jet. 

The Z-sources, which tend to have much longer orbital 
periods (and hence larger orbital separations) than atolls 
and MSXPs, dominate the highest X-ray luminosities. These 
sources spend most of their time in a soft X-ray state and 
in fact have radio luminosities lower than predicted by the 
NSXB hard state radio-X-ra y relation if the ra dio emission 
originates in the jet (Migliari fc Fenderl |2006| ). We would 
therefore expect the OIR emission from the jets in Z-sources 
to also be lower than the model in Fig. 3 unless the radio- 
OIR jet spectrum is in v erted (positive) . In fact from Fig. 
3 of iMigliari fc Fended (|2006l ') we see that at Lx ~ 10^* 
erg s~^, Lradio ~ 10^° erg s~^ and so from Fig. 3 here, the 
radio-OIR spectral index for Z-sourcefl is a ~ 0.2. With 
this information alone, the jets can only dominate the OIR 
of Z-sources if the radio-OIR spectrum is inverted. 



3.2 OIR Spectral Index — X-ray Relations 

Since the dominating OIR emission processes in NSXBs are 
unclear (at least at high Lx) from the OIR-X-ray relations, 
we now turn to the OIR spectra. We see from the lower 
panel of Fig. 3 that there is a relation between Lx and qnir 
for MSXPs and atolls: the NIR spectrum becomes redder at 
higher luminosities. Quantitatively, onir becomes negative 
when Lx > 10^® erg s~^ and cvnir is positive for all data 
below Lx ~ 10^^ erg s~^. This is opposite to the behaviour 
of BHXBs, where there is evidence in some systems for the 
OIR spectru m in th e hard state to become redder at low 
luminosities (iPapeii H. If the origin of the emission is the 
disc blackbody, we would expect a bluer (hotter) spectrum 
at higher luminosities. 

The only process expected to produce an OIR spec- 
trum of index q < at high luminosities in these systems 
is optically thin synchrotron. It is therefore intriguing that 
Qnir < for atolls and MSXPs when Lx > 10^^ ergs"^ 
s~^; the X-ray luminosity range in which the jet could play 
a role (Fig. 3, upper panel). Since there are just five NIR 
data points in the lower panel of Fig. 3 below Lx ~ W^"' erg 



^ At Lx ~ 10'^* erg s"'^ we can calculate the spectral index 
between radio and OIR for Z-sources since we know the radio 
and OIR luminosities. 
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Figure 2. Quasi-simultaneous a (where a is the orbital sepa- 
ration) versus OIR monochromatic luminosity (upper panel) and 
optical and NIR spectral index (lower panel) for NSXBs. OPT 
and NIR refer to the optical and NIR wavebands, respectively 
(also applies to Fig. 3). In the upper panel, the model for opti- 
cal emi ssion from X-ray r eprocessing (solid line) is derived from 
vP&M jvan Paradiis fc N lcClintock 1994) and normalised to the 
optical data of Paper I. The dotted lines show the expected re- 
lation for X-ray reprocessing from NIR data, assuming a NIR- 
optical spectral index of a = 0.5 (upper dotted line) and a = 2.0 
(lower dotted line). In the lower panel, the Z-sources have dif- 
ferent symbols to the atolls and MSXPs (also applies to Fig. 3). 
Data were not included for sources with unknown orbital periods 
or companion masses (Table 1 , and Table 2 of Paper I) . Data for 
Cir X— 1 are denoted by open circles; the orbital separation may 
not be accurate for this source (see Section 2.3). 



s~^, we perform a Kolmogorov-Smirnov (K-S) test to quan- 
tify the significance of the apparent ONiR-ix relation for the 
atolls/MSXPs. We use the 'Numerical recipes in FORTRAN' 
ijPress et al.lll99d ) routine 'kstwo' which is the K-S test for 
two data sets, to determine if the values of qnir differ signif- 
icantly below and above Lx = 10''^ erg s~^. The maximum 
difference between the cumulative distributions is D = 0.91 
with a corresponding probability of P = 5.0 x lO"**; i.e. the 
probability that the NIR spectral index of the data below 
Lx ~ 10^^ erg s~^ belongs to the same population as the 
data above Lx ~ 10"^^ erg s~^ is 0.05 percent. 

In addition, we have carried out a Spearman's Rank 
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Figure 3. Quasi-simultaneous X-ray luminosity versus OIR 
monochromatic luminosity (upper panel) and optical and NIR 
spectral index (lower panel) for NSXBs. In the upper panel, the 
jet model for NSXBs assumes a flat spectrum from radio to OIR 
and adopts the relation Lj-adio ^ 

L]^*' (MfcF MMighari fc Fended 
I2OO6I) . The models for optical and NIR emission from a viscously 
heated disc are also shown, normalised to the optical and NIR 
data, respectively. 



correlation test on the ONiR-Zog Lx and aopT-log Lx data 
of atolls/MSXPs (Table 4). We find an anti-correlation be- 
tween CfNiR and log Lx at the 3.8ct confidence level, support- 
ing the above K-S test results. We do not find a significant 
relation between qopt and log Lx except when we impose 
a somewhat arbitrary X-ray luminosity cut: for data above 
Lx > 10^^ erg s~^ there is an anti-correlation at the 3. So" 
confidence level. This again could be due to the jet contri- 
bution dominating at these highest luminosities. The confi- 
dence of this result should be taken with caution as it could 
be dominated by the group of data with the highest qqpt 
values which happen to lie just above the Lx ~ 10^^ erg s~^ 
cut. 

The optical spectral index qqpt is generally positive 
for atolls and MSXPs (Fig. 3 lower panel), but decreases 
at Lx > 10^^ erg s~^. We would expect aoPT to become 
negative at a higher X-ray luminosity than q:nir if the origin 
of the redder emission component is the jets. As the X-ray 
luminosity is increased, the optically thin synchrotron jet 
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Table 4. Results of the Spearman's Rank correlation between a 
and log Lx for atolls/MSXPs. 



OIR data 


Range in Lx 


Correlation 


Significance 


used 


(erg s-i) 


coefficient 




"NIR 


all 


Ts = -0.63 


3.8a 


aoPT 


all 


rs = -0.14 


O.So- 


aopT 


Lx > 10^'^ 


rs = -0.77 


3. So- 



component will dominate over X-ray reprocessing in the NIR 
bands before the optical as the jet component has a negative 
spectral index. At Lx > 10"^^ erg s~^, oopt is negative 
in 63 percent of the optical data of atolls/MSXPs. Below 
Lx = 10''^ erg s~^ this is 7 percent. These few data in the 
latter group with aoPT < are at low Lx (< 10^* erg s~^) 
and may be du e to cooler accretion discs (possibly like the 
BHXBs; IPaped B). However, the mean uncertainty in the 
values of a are fairly large, so we can only make conclusions 
from general trends and not individual data points. 

Almost all of the spectra of the Z-sources are blue 
(a > 0; Fig. 3 lower panel). This supports the suggestion 
in Section 3.1 that X-ray reprocessing dominates the OIR of 
the Z-sources due to their larger accretion discs and lower 
radio jet luminosities (at a given Lx). The Z-sources cannot 
be dominated by optically thin synchrotron emission as this 
requires a < 0. Since the radio-OIR spectrum of Z-sources 
is a « 0.2 (Section 3.1), the optically thick part of the jet 
spectrum could dominate the OIR if aoiR ~ 0.2 also (al- 
though this requires the optically thick/thin break to be in 
the optical regime or bluer). q:opt,nir > 0.2 is observed for 
most of the data from the Z-sources (Fig. 3 lower panel), 
which implies this is not the case, however we cannot rule 
out an optically thick jet origin to a few of the OIR Z-source 
data (those with lower spectral index). 

The lower panel of Fig. 2 shows how the optical and 

1/2 

NIR spectral index changes with L-^ a. The Z-source data 

1/2 

in this panel have the largest values of L-^ a. The expected 
level of OIR emission from reprocessing in the disc in these 
systems is larger than that of the atolls/MSXPs due to their 
larger accretion discs and higher X-ray luminosities, further 
supporting the X-ray reprocessing scenario. 

3.3 The OIR SEDs 

The SEDs of individual sources in Fig. 4 allow us to see how 
the luminosity-dependent optical and NIR spectral indices 
differ between sources. We group the panels in Fig. 4 into 
the different types of NSXB: MSXPs in the upper panels, 
then atolls, then Z-sources in the lower panels. 

One can visually see a NIR excess {a < 0) joining a 
blue (a > 0) optical spectrum in the SEDs of four of the five 
MSXPs (the 'V'-shape). This is indicative of two separate 
spectral compo nents and is observed i n at least three hard 
state BHXBs (jCorbel fc Fended l2002l : iHoman et all l2005l : 
IPapeil ll). where it is interpreted as the optically thin jet 
spectrum meeting the thermal spectrum of the accretion 
disc. The NIR excess disappears at low luminosities in a 
total of four outbursts of three MSXPs; no NIR excess is 
seen below L^^oir < lO'^** erg s~^ Hz~^ in any MSXP or in 
fact any NSXB. 



Of the five atolls in Fig. 4, one (Aql X-1; the atoll with 
the most data and highest luminosity) has a negative spec- 
tral index; more-so at high luminosity. There is little NIR 
coverage of atolls in the literature (at least at high luminosi- 
ties) and from the data we have, the OIR SEDs are positive 
in all sources except Aql X-1, consistent with the disc dom- 
inating the OIR. The optical spectral index for Aql X-1 is 
very negative at the highest luminosities, suggesting all of 
the OIR was dominated by the jets during the peak of the 
bright 1978 outburst of this source. The SED of 4U 0614-f09 
flattens (becomes redder) in the NIR; this i s known to be 
wher e the disc spectrum meets that of the jet (jMigliari et al.l 
l2006l ). 

There is no evidence for a jet {a < 0) component con- 
tributing to the OIR spectra of any of the seven Z-sources 
in Fig. 4. As was shown above, the jets should be OIR- 
bright in Z-sources but because of their larger discs. X-ray 
reprocessin g dominates. For th e new transient Z-source XTE 
J1701-462 (|Homan et alll2007l ). the optical extinction is not 
well constrained and we show three SEDs in Fig. 4 of the 
same data, using different values of Av to illustrate the pos- 
sible range in spectral indices (we have not included any 
data from XTE J1701-462 in the lower panels of Figs. 2 or 
3 because of these uncertainties in a). 

3.4 The Full Picture 

From the information gathered in Sections 3.1-3.3, the pic- 
ture is now emerging of the dominating OIR emission pro- 
cesses for the different types of NSXB at different luminosi- 
ties. For NSXBs at low luminosities (Lx < 10^® erg s~^), 
the correlations and spectra are consistent with X-ray re- 
processing dominating the OIR. The data are inconsistent 
with emission from the jets at these luminosities. 

We can rule out a viscously heated disc origin to the 
NIR data of all NSXBs since that requires Lnir oc ^x'' 
which is not observed (Fig. 3). Similarly the viscous disc 
cannot dominate the optical light of NSXBs at Lx > lO^'^'^ 
erg s~^ but could below this luminosity as Lqpt oc Lx°'^ is 
required and observed. The jets dominate the NIR and op- 
tical light of atolls/MSXPs above Lx ^ lO^'' and Lx ^ 10^'' 
erg s~^ respectively, as the spectral indices of the contin- 
uum cannot be explained by thermal emission. The SEDs 
of the transient NSXBs with the best sampling show the 
physical disappearance at low luminosity of a NIR excess 
which is present at high luminosity. These observations are 
also consistent with the jet OIR-X-ray model. 

For Z-sources, the OIR-X-ray relations and OIR spectra 
are consistent with emission from X-ray reprocessing. The 
jets and viscously heated disc can be ruled out in most Z- 
sources, however the optical and NIR spectral indices in a 
few observations are also consistent with an optically thick 
jet which extends from the (measured) radio flux. 

We can make a direct measurement of the jet radio- 
NIR spectral index for the atolls/MSXPs, using the NIR 
data which are dominated by the jets. In the luminosity 
range lO^'' < Lx < 10^'^ erg s"\ the NIR data with a 
negative spectral index is on average 0.70 dex more lumi- 
nous than expected from the jet model (which assumes a 
flat radio-NI R spectrum; we know th e radio luminosity at 
this Lx from iMigliari fc Fender! [20061 ). Some low level con- 
tribution from the disc could only partly explain this ex- 
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Figure 4. Spectral Energy Distributions (SEDs) of 17 NSXBs. The Itey in the top left panel corresponds to the slope of the continuum, 
i.e. the spectral index a (where Li, oc jVISXP, ATOLL and Z-SOURCE refer to miUiseco nd X-ray pulsars, atoll sources and Z-sources, 
respectively. The references for the data are listed in Tables 1 and 2 and Table 2 o f lPaped H. The mean systematic error for the data of 
each source is indicated in each panel (in the top left corner in most cases). 



Optical/IR jets in neutron star X-ray binaries 9 



cess. The corresponding radio-NIR spectral index in that 
range of Lx is a ~ 0.16 (0.70 dex in luminosity divided 
by 4.5 dex in frequency between radio and NIR). The opti- 
cally thick to optically thin br eak frequency in NS XBs is 
thought to be in the mid-IR l|Migliari et al.ll20o3 ). mak- 
ing the optically thick radio-mid-IR jet spectrum more in- 
verted; Qthick > 0.2. The spectral evidence for the existence 
of synchrotron-emitting jets in NSXBs presented in this pa - 
per (and showed for the first time bv iMigliari et al.ll2006l ). 
alon g with radio detections which are sometime s resolved 
fsee lFomalont et aLllioOll : iMighari fc Fendeij|2006l ). have di- 
rect implications for the local conditions an d accretion p ro- 
cesses in NSXBs. For example, according to lMassJ l|2006l ). a 
jet can only occur if the magnetic field intensity and mass 
accretion rate are constrained thus: 

/ B m Y'" 

"■^HTo^gJ l^lO-SA/oyr-i ) 

From Table 1 of lMassil (|2006l '). MSXPs only satisfy this 
condition if we adopt the largest mass accretion rates and 
smallest magnetic field strengths measured. Empirically we 
do see the jet only at high luminosities and hence high mass 
accretion rates. The atolls and Z-sources more easily satisfy 
the condition for jet production, except for the highest mag- 
netic field strengths measured in Z-sources. Since we know 
jets exist in all three flavours of NSXB at least at high lu- 
minosities, perhaps the measured mass accretion rates and 
magnetic field strengths are under- and over-estimated re- 
spectively, in some cases. 



4 CONCLUSIONS 

We have shown that the dominating OIR emission processes 
in NSXBs vary with X-ray luminosity and between source 
types (atoUs/MSXPs and Z-sources). However, a clearer pic- 
ture seems to em erge th an the dominating emission pro- 
cesses of BHXBs (jPaperl 13) . Models predict that X-ray re- 
processing in the accretion disc should dominate the OIR at 
low luminosities and the jets, if present, should dominate at 
high luminosities for NSXBs with relatively small accretion 
discs. 

From the observed spectral index of the OIR continuum 
and from OIR-X-ray relations we show that this is the case 
in atolls and MSXPs: the jets dominate the NIR and optical 
emission above Lx ~ 10'^'' erg s~^ and Lx ~ 10'^^ erg s~^, re- 
spectively. Below these luminosities X-ray reprocessing dom- 
inates, although we cannot rule out a viscously heated disc 
origin to the optical emission. We have shown that the radio- 
NIR spectral index of the jets in the atoUs/MSXPs is slightly 
inverted: a ~ 0.16, at least at high X-ray luminosities (10''^- 
10''^ erg s~^). In the Z-sources, which have larger discs, we 
find that X-ray reprocessing is responsible for all OIR emis- 
sion, and the radio-OIR jet spectrum has to be a < 0.2 
(otherwise the jet spectrum would dominate over the disc). 
However, the optically thick part of the jet spectrum could 
dominate in a few cases. 

Evidence for the existence of NSXB jets is mounting 
both fr om radio detections (which are sometimes resolved; 
see e.g. iFomalont et al.|[200ll ) and spectrally from higher fre- 
quency observations (Migliari et al. 2006; this paper). The 



power carried in NSXB jets is sensitive to the position of the 
break between optically thick and thin emission in its spec- 
trum, and is curr ently a topic of deba te. This is constrained 
in 4U 0614-1-09 (|Mighari et al.ll200d ). where it is likely to 
lie in the mid-infrared. Mid-infrared and NIR photometry 
and polarimetry could shed light on this; optically thin syn- 
chrotron emission is expected to be highly polarised if the 
magnetic field is ordered. 
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